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The  ONION  method  is  a  reflection-coefficient  measurement  technique  designed  for  use 
on  data  acquired  from  thick  underwater  acoustic  panels  in  the  frequency  range 
1-10  kHz,  but  may  be  used  to  frequencies  as  high  as  25  kHz.  The  method  extrapolates 
transient  ref lected-wave  data  using  least-squares  fitting  to  a  multiple-layer  panel 
model.  A  description  of  the  method,  as  it  applies  to  data  acquired  using  a  normal¬ 
ly  incident  interrogating  wave,  is  provided  in  J".U.  PiqUeEt®^  *The  ONION  method: 

A  reflection  coefficient  measurement  technique  for  thick  underwater  acoustic 
panels,  1  J -- Acoust-.— Socr-Amr' 85, —1029- 1040  (1989)^/  The  present  article  describes 
research  directed  toward  generalizing  the  method  to  allow  for  an  offnormal  incidence 
interrogating  wave.  Successful  applications  of  the  generalized  method  to  data  ac¬ 
quired  under  offnormal  incidence  are  described.  The  generalized  method  was  applied 
to  measurements  that  were  made  on  two  different  sample  panels.  The  measurements 
involved  one  test  temperature,  two  hydrostatic  test  pressures,  and  three  test:  fre¬ 
quencies  .  '  i 
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The  ONION  method  is  a  reflection-coefficient  measurement  technique  designed  for  use  on 
data  acquired  from  thick  underwater  acoustic  panels  in  the  frequency  range  1-10  kHz,  but 
may  be  used  to  frequencies  as  high  as  25  kHz.  The  method  extrapolates  transient  reflected- 
wave  data  using  least-squares  fitting  to  a  multiple-layer  panel  model.  A  description  of  the 
method,  as  it  applies  to  data  acquired  using  a  normally  incident  interrogating  wave,  is  provided 
in  J.  C.  Piquette,  "The  ONION  method:  A  reflection  coefficient  measurement  technique  for 
thick  underwater  acoustic  panels,”  J.  Acoust.  Soc.  Am.  85,  1029-1040  ( 1989).  The  present 
article  describes  research  directed  toward  generalizing  the  method  to  allow  for  an  offnormal 
incidence  interrogating  wave.  Successful  applications  of  the  generalized  method  to  data 
acquired  under  offnormal  incidence  are  described  The  generalized  method  was  applied  to 
measurements  that  were  made  on  two  different  sample  panels.  The  measurements  involved  one 
test  temperature,  two  hydrostatic  test  pressures,  and  three  test  frequencies. 

PACS  numbers:  43.20.Fn,  43.20.Px,  43.30.Sf,  43.60  Gk 


INTRODUCTION 

Panel  measurements  are  a  standard  technique  w  hereby 
the  effectiveness  of  a  coating  material  at  reducing  unwanted 
echoes  is  determined.  The  conventional  panel-measurement 
method1 :  involves  the  use  of  a  sample  panel  whose  lateral 
dimensions  are  large  compared  to  a  wavelength  of  the  inter¬ 
rogating  wave  in  the  surrounding  fluid  medium.  However, 
due  to  the  difficulty  and  cost  of  fabricating  large  samples, 
and  due  to  the  limited  size  of  test  facilities,  available  samples 
often  have  lateral  dimensions  that  are  less  than  one  such 
wavelength  for  frequency  ranges  of  interest. 

In  order  to  avoid  the  interfering  influence  of  the  diffract¬ 
ed  waves  originating  at  the  sample  edges,  it  is  often  necessary 
to  operate  a  panel  test  in  the  pulsed  mode,  and  to  utilize 
portions  of  the  experimentally  measured  reflected-wave 
pulses  that  have  not  achieved  steady  state.  One  technique  for 
treating  such  nonsteady-state  reflected  wave  signals,  which 
is  applicable  when  the  observed  transient  behaves  as  if 
caused  by  a  lumped-parameter  system,  is  the  Prony'  meth¬ 
od.  Such  behavior  arises  primarily  in  panels  of  small  overall 
thickness,  so  that  the  majority  of  the  reflected-wave  tran¬ 
sient  is  caused  by  the  turnon  transient  of  the  source  of  the 
interrogating  wave. 

As  panel  thickness  is  increased,  the  sample  behaves  in¬ 
creasingly  as  a  distributed-parameter  system,  so  that  the 
Prony  method  is  no  longer  applicable  to  the  reflected  wave¬ 
form.  The  ONION  method4  5  is  ideally  suited  to  samples 
that  behave  as  a  distributed-parameter  system;  i.e.,  samples 
whose  reflected  wave  contains  significant  transients  that  are 
associated  with  the  round-trip  travel  times  of  waves  in  the 
panel  subla'"'":  This  method  is  a  reflection-coefficient  mea¬ 
surement  technique  that  is  designed  for  use  on  thick  under¬ 


water  acoustic  panels  in  the  frequency  range  1-10  kHz,  but 
may  be  used  to  frequencies  as  high  as  25  kHz.  The  approach 
is  based  on  least-squares  fitting  of  a  multiple-layer  panel 
model  to  the  measured  reflected-pulse  waveform.  For  a  com¬ 
plete  description  of  die  method,  see  Ref.  4. 

The  initial  development  of  the  ONION  method  was  re¬ 
stricted  to  the  analysis  of  reflected-wave  data  acquired  by 
probing  the  sample  with  an  interrogating  wave  that  arrives 
at  normal  incidence  to  the  panel.  However,  since  the  reflec¬ 
tion  coefficient  is  a  function  of  incidence  angle,  it  is  of  inter¬ 
est  to  determine  sample  behavior  as  a  function  of  measure¬ 
ment  angle.  The  present  article  examines  a  generalization  of 
the  ONION  method  to  allow  for  offnormal  incidence  angles. 

Section  I  presents  a  synopsis  of  the  approach.  A  discus¬ 
sion  of  geometry,  edge  waves,  and  certain  ramifications  of 
the  theoretical  panel  model  is  presented  in  Sec.  II.  This  sec¬ 
tion  also  describes  how  angular  interpolations  and  extrapo¬ 
lations  are  achieved.  A  description  of  experimental  measure¬ 
ments  made  to  investigate  the  effectiveness  of  the  method  is 
given  in  Sec.  III.  Section  IV  gives  a  discussion  of  the  meaning 
of  the  measurements  obtained,  and  also  presents  a  descrip¬ 
tion  of  some  potential  influences  of  certain  experimental 
aspects  that  represent  departures  from  the  ideal  conditions 
assumed  by  theory.  A  summary  and  the  conclusions  are  giv¬ 
en  in  Sec.  V. 

Throughout  this  paper,  it  is  assumed  that  the  reader  is 
thoroughly  familiar  with  the  ONION  method  as  applied  to 
the  normal-incidence  case.  It  is  also  assumed  that  the  reader 
understands  the  structure  and  behavior  of  an  underwater 
acoustic  panel,  and  is  familiar  with  the  conventional  panel- 
measurement  configuration.  The  reader  who  is  not  confident 
in  his  knowledge  of  these  subjects  is  directed  to  Refs.  1 ,  2, 4, 
and  5. 
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I.  BRIEF  OVERVIEW  OF  THE  GENERALIZED  ONION 
METHOD 

Figure  1  is  a  block  diagram  that  summarizes  the  ap¬ 
proach  used  in  the  generalized  ONION  method.  The  algo¬ 
rithm  involves  three  phases.  In  phase  1.  incident  and  reflect¬ 
ed  pulsed  waveforms  experimentally  measured  at  normal 
incidence  are  used  in  a  normal-incidence  ONION-method 
calculation'  to  obtain  starting  model  parameters.  In  phase  2. 
a  further  least-squares  adjustment  of  the  model  parameters 
is  performed.  During  this  phase,  a  simultaneous  least- 
squares  minimization  of  the  mean-squared  error  between  the 
model  and  the  data  is  performed  using  all  incident  and  re¬ 
flected  pulsed  waveforms  that  hav  e  been  measured  at  each  of 
the  incidence  angles  of  interest.  The  maximum  permissible 
incidence  angle  is  restricted  by  the  requirement  to  avoid  the 
interfering  influence  of  the  sample  edge-difTracted  wave. 

( This  will  be  described  further  in  Sec.  II  B. )  Since  a  restrict¬ 
ed  range  of  incidence  angles  is  considered,  it  is  assumed  that 
a  //w/V/-layer  panel  model  may  be  used  in  the  fitting  process. 
(Justifications  for  this  assumption  are  presented  in  Sec. 
IIC.)  Phase  2  is  iterated  several  times  to  achieve  model  pa¬ 
rameters  that  are  most  consistent  (in  a  least-squares  sense) 
with  the  experimental  data. 

Once  the  "best-fit"  parameter  values  have  been  deter¬ 
mined  by  phase  2,  these  parameter  values  are  substituted 
into  a  vo//<r/-layer  panel  model,  during  phase  3.  to  deduce 
reflection  coefficients  throughout  the  range  of  incidence  an¬ 
gles  from  0  to  89  deg.  (The  90-deg  direction  is  omitted  to 
avoid  a  numerical  singularity  in  the  available  software  that 
performs  the  calculation. )  In  using  the  solid-layer  model,  a 
priori  values  of  the  shear  properties  of  the  layers  are  as¬ 
sumed.  In  most  cases  of  measurement  interest,  however,  it  is 
assumed  that  the  materials  of  the  panel  layers  are  characteri- 
zable  by  negligible  shear;  it  is  assumed  that  only  the  steel 
backing  plate  contributes  a  significant  shear  effect.  The 
phase  3  portion  of  the  algorithm,  which  implements  the  sol¬ 
id-layer  panel  model  calculation  just  described,  computes 
reflection  coefficients  that  interpolate  between  measurement 
angles  and  that  extrapolate  beyond  the  greatest  available 
measurement  angle. 

II.  CALCULATIONAL  CONSIDERATIONS 
A.  Geometry 

Figure  2  depicts  the  measurement  configuration  used 
for  acquiring  refiected-wave  data  in  an  offnormal-incidence 
panel  measurement.  In  this  figure,  d ,  is  the  distance  separat¬ 
ing  the  source  and  the  rotator  shaft  that  supports  the  test 
panel.  The  quantity  </,,  represents  the  offset  distance  of  the 
d-mrtor  hydrophone  from  the  first  layer,  i.e.,  the  layer  clos¬ 
est  to  the  hydrophone.  The  hydrophone  is  rigged  so  that  it 
corotates  with  the  test  panel;  i.e..  it  remains  centered  with 
respect  to  the  panel’s  edges.  The  quantity  t  is  the  overall 
sample  thickness.  The  angled,  which  is  depicted  in  the  figure 
as  the  angular  location  of  the  hydrophone  with  respect  to  the 
acoustic  axis,  also  repres^n*'  the  rotator  shaft  anele  that  is 
adjusted  during  measurement.  The  angle  0  represents  the 
incidence  angle  for  specularly  reflected  waves. 

It  is  important  to  note  that  the  angles  6  and  4>  are  gener¬ 
ally  unequal.  A  trigonometric  calculation  shows  that 
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FIG.  !.  Block  diagram  of  the  generalized  ONION-melhod  algorithm.  The 
algorithm  involves  three  phases.  In  phase  1.  the  experimentally  measured 
incident  pulsed  waveform  at  normal  incidence p\  ( ; ) .  and  the  resulting  ex¬ 
perimentally  measured  reflected  pulsed  waveform  p'"'U).  are  used  in  a  nor¬ 
mal-incidence  ONION-method  calculation  to  obtain  starting  model  pa¬ 
rameters.  In  phase  2.  the  model  parameters  obtained  from  phase  1  are 
iteratively  improved  using  a  nonlinear  least-squares  fitting  procedure  that 
simultaneously  fits  an  oflhormal  incidence  theoretical  panel  model  to  data 
acquired  at  all  measurement  angles.  A  fluid- layer  panel  model  is  used  dur¬ 
ing  this  phase.  Finally,  in  phase  3.  the  best-fit  model  parameters  deduced  by 
phase  2  are  used  together  with  assumed  shear  properties  for  the  layers  in  a 
calculeuor  based  on  a  to/fo'-layer  panel  model  to  obtain  interpolated  and 
extrap  e  .  '’flection  coefficients  asafunction  of  incidence  angle  8.  Here. 
p\"'(i )  ‘s  the  experimentally-measured  incident  pulsed  time  wave¬ 

form  at  ih  incidence  angle  (n  =  0  represents  normal  incidence); 
pl"'U)  represents  the  experimentally  measured  reflected  pulsed  time  wave¬ 
form  at  the  nth  incidence  angle;  andp)’"(()  is  the  computed  reflected  pulsed 
time  waveform  at  the  nth  incidence  angle,  based  on  a  multiple  fluid-layer 
panel  model. 


9  =  tan 


t 


d,  sin  4> 


;)• 


(i) 


\ds  cos  <t>  -  /  +dh/ 

For  a  panel  nic.-urtnient  as  performed  in  the  Anechoic 
Tank  Facility  (ATF)  of  the  Underwater  Sound  Reference 
Detachment  of  the  Naval  Research  Laboratory  (NRL- 
USRD)  in  OrU  "do.  Florida,  a  conventional  couiiguratton 
uses  values  of  d,  =  170  cm  and  dh  =  34  cm.  For  a  panel  of 
thickness  t  =  15  cm  and  a  rotator  shaft  measurement  angle 
of  <j>  =  45°,  Eq.  ( 1 )  yields  an  incidence  angle  of  40.8°; 
thus  a  considerable  angular  error  can  be  made  if  tfand  <t>  are 
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FIG.  2.  Diagram  depicting  the  geometry  of  an  offnormal-incidence  panel 
measurement  configuration.  The  figure  presents  a  top  view;  i.e.,  from  di¬ 
rectly  above  the  test  panel.  Here,  rf,  represents  the  source-to-rotator  shaft 
separation  distance;  dh  represents  the  hydrophone  offset  distance  from  the 
first  panel  layer  (i.e.,  the  panel  layer  closest  to  the  hydrophone);  t  repre¬ 
sents  the  total  panel  thickness;  6  represents  the  rotator  shaft  angle;  and  6 
represents  the  incidence  angle  for  specular  reflection. 


assumed  to  be  identical.  It  is  interesting  to  note,  however, 
that  if  the  hydrophone  offset  is  equal  to  the  panel  thickness, 
i.e..  if  dh  =  r,  then  Eq.  ( 1 )  reduces  to  6  —  d>  for  all  measure¬ 
ment  angles. 

T  ypical  values  of  the  measurement  variables  of  Fig.  2  as 
used  in  ONION-method  measurements  are  d%  =  200  cm 
and  dh  =  5  cm.  For  a  panel  of  t  —  15  cm  overall  thickness 
and  a  measurement  angle  of  <b  =  45°,  Eq.  ( 1 )  leads  to  the 
value  0  —  47. 1°.  Note  that  ^exceeds  for  the  typical  configu¬ 
ration  used  in  ONION-method  measurements,  while  9  is 
less  than  <t>  for  the  typical  configuration  used  in  conventional 
measurements.  This  difference  in  behavior  is  primarily  due 
to  the  different  hydrophone  offsets  dh  used  in  each  configu¬ 
ration. 

It  is  also  important  to  realize  that  for  panels  containing 
layers  composed  of  deformable  materials,  the  overall  panel 
thickness  t  will  vary  as  a  function  of  hydrostatic  test  pres¬ 
sure.  (The  hydrophone  offset^  will  also  change;  this  will  be 
discussed  further  in  the  following  subsection. )  Thus  Eq.  ( 1 ) 
will  yield  a  different  incidence  angle  9  for  the  same  rotator 
shaft  angle  <b  at  different  hydrostatic  test  pressures. 

B.  Edge  wave®  and  data  windows 

Since  the  theoretical  model  used  to  deduce  the  reflected 
pulsed  waveform  in  the  ONION  method  neglects  wave  con¬ 
tributions  arising  from  the  sample  edges,  it  is  important  to 
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FIG.  3.  Diagram  depicting  acoustic  ray  paths  for  the  specularly  reflected 
wave  and  the  edge-diffracted  wave  arising  from  the  steel  backing  plate.  The 
figure  presents  a  top  view,  i.e.,  from  directly  above  the  test  panel.  Variables 
are  the  same  as  those  used  in  Fig.  2.  The  additional  variable  w  represents  the 
width  of  the  sample  panel  (assumed  to  have  a  square  cross  section). 


choose  data  “windows”  that  admit  into  the  analysis  only 
data  that  are  not  significantly  corrupted  by  panel-edge  dif¬ 
fraction.  Unfortunately,  due  to  the  relatively  low  sound 
speeds  and  relatively  large  thicknesses  of  samples  of  mea¬ 
surement  interest,  it  is  usually  not  possible  to  entirely  ex¬ 
clude  contributions  from  panel-edge  waves.  However,  this 
problem  is  not  severe  due  to  the  fact  that  the  initial  sample 
layers  (i.e.,  thoseclosest  to  the  hydrophone)  typically  havea 
good  acoustic  “match"  to  the  surrounding  water  medium, 
and  the  interrogating  wave  has  a  turnon  transient;  hence, 
initial  edge-wave  contributions  tend  to  be  of  rather  low  am¬ 
plitude.  On  the  other  hand,  the  edge  wave  due  to  the  steel 
backing  plate  that  is  typically  affixed  to  test  panels  can  be 
rather  large  in  amplitude;  thus  the  data  window  is  selected  to 
avoid  edge-wave  contributions  associated  with  the  backing 
plate. 

In  view  of  the  fact  that  the  materials  used  to  fabricate 
sample  panels  are  typically  characterized  by  sound  speeds 
less  than  that  of  the  surrounding  water  medium,  it  is  usually 
the  acoustic  path  through  t'ie  water  medium  that  represents 
the  least  time-of-flight  between  any  two  points  of  interest  in  a 
panel  measurement.  A  typical  source-to-backing  plate-to- 
detector  path,  and  the  specularly  reflected  path,  are  depicted 
in  Fig.  3. 


las 

r 


Jean  C.  Piquette:  Offnormal  reflection  coefficient 


In  terms  of  the  variables  of  Fig.  3.  we  can  deduce  expres¬ 
sions  for  the  path  length  of  the  specularly  reflected  wave 
and  the  path  length  of  the  edge-diffracted  wave 
arising  from  the  backing  plate.  The  expressions  are 

=  [  (^/,  +  d,  cos  <t>  —  /)-  +  (ds  sin  <£):] (2) 

and 

=  {[  (tt/2)cos  <b\-  +  [/f,  —  f  tt/2 )sin  ]  - } 1  2  +  / 

-t-  [  ( w/2 )  ~  4-  d  ~h  ] 1  ",  ( 3 ) 

where  w  represents  the  sample  width  ( the  sample  is  assumed 
to  have  a  square  cross  section).  In  applying  Eqs.  (2)  and 
( 3),  it  should  be  realized  that  when  hydrostatic  test  pressure 
is  increased,  overall  panel  thickness”  t  decreases,  while  the 
hydrophone  offset  undergoes  a  corresponding  increase. 
This  is  so  because  the  hydrophone  is  rigged  to  maintain  a 
fixed  distance  from  the  rotator  shaft,  but  the  quantity  dh  is 
measured  with  respect  to  the  front  panel  face.  This  panel  face 
recedes  from  the  hydrophone  as  hydrostatic  test  pressure  is 
increased.  The  two  effects  (i.e.,  increasing  d„  and  decreasing 
/)  tend  to  decrease  the  acoustic  path  difference  6^,  -  5  . 

Hence,  if  a  single  data  window'  width  is  to  be  used  at  a  given 
measurement  angle  for  all  hydrostatic  test  pressures,  a  con¬ 
servative  approach  is  to  use  the  window'  width  that  elimi¬ 
nates  the  edge-wave  contribution  at  the  greatest  hydrostatic 
test  pressure  of  interest. 

For  example,  imagine  a  sample  panel  whose  overall 
thickness  /  is  15  cm  at  atmospheric  pressure,  but  which  de¬ 
creases  to  12.5  cm  at  the  greatest  hydrostatic  test  pressure  of 
interest.  (Such  a  compression  is  realistic  for  the  panels  of 
measurement  interest.)  Imagine  also  that  the  rotator  shaft 
angle  of  interest  is  6  =  20°.  It  is  convenient  in  the  calcula¬ 
tions  that  follow  to  imagine  that  the  incident  and  reflected 
wave  data  are  acquired  in  digital  form,  as  required  by  the 
ONION-method  algorithm.  Let  us  take  the  digital  data  ac¬ 
quisition  rate  to  be  4  MHz,  and  calculate  the  corresponding 
number  of  reflected-wave  data  points  in  the  data  window 
that  avoids  contributions  from  the  edge-diffracted  wave  aris¬ 
ing  from  the  steel  backing  plate.  This  number  of  points  is 
calculated  from  the  formula 


where  c„  is  the  speed  of  sound  in  the  surrounding  water 
medium.  Assuming  a  hydrophone  offset  of  dh  —  5  cm  at  at¬ 
mospheric  pressure,  a  water  sound  speed  of  c0  =  1.5  X  10' 
cm/s,  and  a  4-MHz  data  measurement  rate,  Eq.  (4)  yields  a 
data  window  width  of  1435  points  at  <f>  =  20°,  assuming  a 
sample  width  of  w  =  76cm  (30  in.).  However,  recalculating 
the  data  window  width  using  the  !2.5-cm  thickness  assumed 
for  the  highei  hydrostatic  test  pressure  results  in  a  window 
width  of  1 250  points.  ( In  this  calculation,  a  hydrophone  off¬ 
set  of  dh  =  7.5  cm  must  be  used;  i.e,  the  5-cm  offset  given  for 
atmospheric  pressure  must  be  increased  by  an  amount  equal 
to  the  corresponding  decrease  in  sample  thickness,  which  is 
15  —  12.5  =  2.5  cm  in  tins  example.)  Thus  to  avoid7  the 
edge  wave  from  the  backing  plate  at  all  hydrostatic  test  pres¬ 
sures  of  interest  requires  a  window  width  of  no  greater  than 
1 250  points  at  the  rotator  shaft  angle  <t>  =  20°. 


C.  Model  considerations 

The  theoretical  model  used  in  the  normal- incidence  im¬ 
plementation  of  the  ON  ION  method4  s  treats  each  panel  lay¬ 
er  as  a  fluid;  i.e.,  it  is  assumed  tha*  the  shear  modulus  of  each 
layer  is  negligible.  This  simplification  is  readily  justifiable  in 
the  normal-incidence  case.  Panel  layers  in  samples  of  inter¬ 
est  are  generally  fabricated  in  a  manner  that  creates  in  each 
layer  an  array  of  air-filled  macrovoids  in  a  rubber  matrix 
material.  Both  air  and  rubber  are  characterized  by  quite  low 
values  of  the  shear  modulus.  Of  course,  the  test  panel  is  typi¬ 
cally  affixed  to  a  support  plate  fabricated  from  steel,  a  mate¬ 
rial  having  a  very  significant  shear  modulus.  However,  it 
should  be  noted  that  no  shear  waves  are  excited  in  a  solid  of 
infinite  lateral  extent  that  is  stimulated  by  a  u.'nnally  inci¬ 
dent  interrogating  wave.  The  influence  of  the  edges  of  a  finite 
sample  can  be  reduced  by  acquiring  the  measurements  in  the 
pulsed  mode,  and  gating  out  edge  effects.  Hence,  even  the 
steel  backing  plate  can  accurately  be  modeled  as  a  fluid  in  the 
normal-incidence  case. 

As  the  interrogating  wave  incidence  angle  is  allowed  to 
deviate  from  the  normal,  however,  shear-wave  production 
becomes  increasingly  significant,  so  that  it  is  possible  that 
the  steel  backing  plate  may  no  longer  be  accurately  treated  as 
a  fluid.  (  Due  to  the  very  low  values  of  the  shear  modulus  of 
the  materials  constituting  the  maerovoided  viscoelastic  lay¬ 
ers.  however,  it  is  assumed  that  shear-wave  production  will 
remain  insignificant  for  these  layers  even  at  large  offnormal 
incidence  angles. ) 

The  influence  of  shear-wave  production  in  a  steel  back¬ 
ing  plate  can  be  seen  by  examining  Table  I.  This  table  pre¬ 
sents  complex  reflection  coefficients  as  a  function  of  inci¬ 
dence  angle  for  a  single  layer  of  steel  and  for  a  single  layer  of 
a  steel-like  fluid,  fin  this  table,  and  the  following  tables,  the 
real  and  imaginary  parts  of  the  complex  reflection  coeffi¬ 
cient  are  presented  using  the  complex-number  convention 
(real,  imaginary).]  The  steel-like  fluid  is  taken  to  have  a 


TABLE  I.  Complex  reflection  coefficients  as  a  function  of  incidence  angle 
for  a  single  layer  of  steel  and  for  a  single  layer  of  a  steel-like  fluid.  The  steel¬ 
like  fluid  is  taken  to  have  a  density  equal  to  that  of  steel  and  to  have  a  longi¬ 
tudinal  sound  speed  equal  to  the  longitudinal  wave  speed  in  steel.  In  both 
cases,  the  sample  is  immersed  in  an  infinite  water  medium.  Frequency  is  20 
kHz.  Sample  thickness  is 0.95  cm  (,|  in.). 


Comple: 

v  reflection 

coefficient  (dimensionless) 

Incidence 

angle  0 

Steel-like 

(degrees) 

Steel 

fluid 

0 

0.906,0.289 

0.906,0.289 

5 

0.905,0.290 

0.905,0.290 

10 

0.904,0.292 

0.903.0.294 

15 

0.900.0.304 

0.899.0.301 

20 

0.895,0.301 

0.894.0.31 1 

25 

0.887.0.312 

0.886.0.325 

50 

0.873.0.327 

0875.0.341 

35 

0.855.0.348 

0.861.0.362 

40 

0.826,0.374 

0.843,0.386 

45 

0.785.0.405 

0.819,0.415 
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negligible  shear  modulus  (corresponding  to  a  shear  wave- 
speed  of  1  m/s)  and  is  taken  to  have  the  same  density  as  steel 
and  to  have  a  longitudinal  sound  speed  equal  to  the  longitu¬ 
dinal  wavespeed  in  steel  In  both  cases,  the  layer  is  assumed 
to  be  immersed  in  an  infinite  water  medium.  Frequency  is  20 
kHz  and  sample  thickness  is  0.95  cm  in  ).  This  thickness 
corresponds  to  the  thickness  of  a  standard  support  plate 
typically  used  in  a  panel  test.  The  20-kIIz.  frequency  consid¬ 
ered  is  close  to  the  maximum  25-kHz.  frequency  for  which 
ONION-method  measurements  are  intended. 

Perhaps  the  most  striking  feature  of  Table  l  is  that  the 
values  presented  for  steel  actually  deviate  very  little  from 
those  presented  for  the  steel  like  fluid,  at  least  for  incidence 
angles  less  than  30°.  This  behavior  is  also  typical  of  that  seen 
for  the  lower  test  frequencies  of  5  and  10  kHz  that  were 
considered,  but  which  are  not  presented  here.  (In  fact,  the 
deviation  is  even  less  for  frequencies  less  than  20  kHz. )  This 
behav  ior  is  probably  caused  by  the  small  layer  thickness  of 
0.95  cm  that  was  used  in  the  calculation.  However,  since  the 
thickness  used  in  the  calculations  that  generated  Table  1  rep¬ 
resents  the  thickness  of  a  standard  support  plate,  and  in  view 
of  the  rather  low  shear  modulus  characteristic  of  the  materi¬ 
als  that  constitute  the  layers  of  test  panels,  the  results  pre¬ 
sented  in  Table  I  are  taken  as  a  justification  of  the  continued 
use  of  a  fluid-layer  model  to  analyze  offnormal  incidence 
measurements  in  the  phase  2  portion  of  the  generalized  algo¬ 
rithm.  at  least  for  incidence  angles  0  less  than  30°.  ( In  fact.  0 
is  restricted  to  be  less  than  or  equal  to  that  incidence  angle 
that  corresponds  to  a  rotator  shaft  angle  of  <]>  =  20°  in  the 
calculations  in  which  a  fluid-layer  model  is  used.  This  angu¬ 
lar  restriction  is  used  to  avoid  the  edge  wave  from  the  steel 
backing  plate. )  The  motivation  for  using  a  fluid-layer  model 
rather  than  the  more  accurate  solid-layer  model  during  the 
phase  2  portion  of  the  method  is  to  avoid  the  significant 
increase  in  computer  processing  time  that  would  result  from 
using  a  solid-layer  model/ 

The  validity  of  the  fluid-layer  model  for  offnormal  inci¬ 
dence  can  be  further  investigated  by  referring  to  the  results 
presented  in  Tables  1 1— IV.  These  tables  present  information 
similar  to  that  presented  in  Table  I,  except  in  these  tables  the 


TABLE  II  Complex  reflection  coefficients  as  a  function  of  incidence  angle 
for  a  simple  three-layc  sample  Frequency  is  5  kHz. 


Complex  reflection 
coefficient  (dimensionless) 


Incidence 
angle  0 
(degrees) 

Solid 

layers 

Fluid 

layers 

0 

O '24,0.218 

0.524,0.218 

5 

0.525,0.211 

0.527,0.211 

10 

0.510,0.200 

0.535,0.189 

15 

0.542,0.169 

0.548,0.152 

20 

0.491,0.218 

0.560,0.100 

25 

0.521,0  161 

0  569,0.033 

.10 

0.521.0. 1 16 

0.570.0.047 

15 

0.513,0.080 

0.558.0.138 

40 

0.496,0.045 

0.528,0.235 

45 

0.470,0.014 

0.477,0.335 

TABLE  III  Complex  reflection  coefficients  as  a  function  of  incidence  angle 
for  a  simple  three-layer  sample.  Frequency  is  10  kHz. 


Complex  reflection 

coefficient 

(dimensionless) 

Incidence 

angle  0 

Solid 

Fluid 

t degrees) 

layers 

layers 

0 

0.582.0.509 

0.582.0.509 

5 

0.587.0.497 

0  591,0.492 

10 

0.603.0.457 

0.616.0.416 

15 

0.627.0.381 

0.648.0.116 

20 

0.641.0.261 

0.666.0.1X2 

25 

0.671.0.126 

0.611.0.029 

30 

0.641.0.057 

0.496.0.270 

25 

0.565,0.218 

0.206.0.458 

40 

0.448.0.191 

0.192.0.457 

45 

0.314.0.494 

0.527.0.206 

sample  considered  consists  of  three  homogeneous  sublayers'1 
instead  of  the  single  steel  layer  considered  in  Table  1.  In  this 
sample,  the  layer  closest  to  the  sound  source  is  composed  of 
polymethylmethacrylate  (PMM)  of2.54-cm  (1-in.)  thick¬ 
ness.  the  second  layer  is  a  water  layer  of  2.54-cm  (1-in.) 
thickness,  and  the  third  layer  is  steel  layer  of0.95-cnt  (cin. ) 
thickness.  The  column  labeled  "solid  layers"  presents  com¬ 
plex  reflection  coefficients  evaluated  by  treating  the  PMM 
and  steel  layers  as  solids.  The  results  presented  in  the  column 
labeled  “fluid  layers"  were  obtained  by  treating  the  PMM 
and  steel  layers  as  having  a  negligible  shear  modulus.  In  ad¬ 
dition.  the  density  of  each  of  these  layers  was  taken  equal  to 
that  in  the  corresponding  solid  layer,  and  the  longitudinal 
sound  speed  in  each  of  these  layers  was  taken  equal  to  the 
longitudinal  wave  speed  of  the  corresponding  solid  layer. 

As  can  be  seen  by  referring  to  Tables  II  and  III,  the 
results  for  the  solid-layer  and  fluid-layer  cases  are  very  simi¬ 
lar  to  each  other  for  the  frequencies  5  and  10  kHz.  at  least  for 
incidence  angles  less  than  or  equal  to  20°.  The  deviation  is 
greater  for  the  20-kHz  case,  as  can  be  seen  by  referring  to 
Table  IV,  and  is  particularly  so  for  9  —  15°.  Note,  however, 


TABLE  IV.  Complex  reflection  coefficients  as  a  function  of  incidence  angle 
for  a  simple  three-layer  sample.  Frequency  is  20  kHz. 


Complex  reflection 
coefficient  (dimensionless) 

Incidence 

angle  6 

Solid 

Fluid 

( degrees ) 

layers 

layers 

0 

0  974.0.164 

0.974.0  164 

5 

0.972,0.172 

0.971.0.182 

10 

0.966,0.203 

0.960.0.236 

15 

0.224,0.970 

0.936,0.319 

20 

0.937,0.305 

0.893,0.425 

25 

0.871,0.459 

0.825,0.546 

30 

0.694,0.697 

0.727,0.671 

35 

0.181,0.970 

0.597,0.788 

40 

0.858.0.512 

0.431,0.887 

45 

0.462,0.846 

0.224.0,954 
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that  this  sample  represents  a  rather  severe  test  of  the  fluid- 
layer  model,  in  view  of  the  significant  shear  modulus  and 
thickness  of  the  PMM  layer.  If  the  shear  modulus  of  the 
PMM  layer  is  permitted  to  decrease  to  an  insignificant  value, 
w  hile  retaining  the  true  shear  modulus  in  the  steel  layer  and 
the  true  longitudinal  wave  speed  for  PMM.  agreement  with 
the  fluid-layer  case  is  substantially  improved.  (This  calcula¬ 
tion  i>-  of  interest  since  the  materials  used  in  panels  of  actual 
measurement  interest  have  a  substantially  lower  shear  mod¬ 
ulus  than  that  of  PMM. )  Table  V  presents  results  for  the  20- 
kH/  case  in  which  the  PMM  layer  is  treated  as  a  fluid  but  the 
steel  layer  ts  treated  as  a  solid.  As  can  be  seen  by  comparing 
Table  V  with  the  fluid-layer  column  of  Table  IV.  agreement 
is  considerably  improved,  even  at  0  —  1 5°.  In  view  of  the  fact 
that  the  situation  considered  in  generating  the  results  dis¬ 
plac'd  in  Table  V’  more  closely  represents  the  situation  of 
actual  measurement  interest  than  the  situation  considered  in 
generating  the  solid-layer  results  of  Table  IV.  we  take  the 
above-mentioned  agreements  as  establishing  the  validity  of 
the  fluid-layer  model,  at  least  for  frequencies  less  than  or 
equal  to  20  kHz  and  for  incidence  angles  corresponding  to 
rotator-shaft  angles  that  are  less  than  or  equal  to  20°.  The 
validity  of  this  assumption  is  also  further  investigated  in  the 
description  of  the  experiments  presented  in  Sec.  III. 

D.  Angular  interpolation  and  extrapolation 

Based  on  t lie  calculations  described  in  the  preceding 
subsection,  and  based  on  the  verifying  measurements  to  be 
described  in  Sec.  III.it  is  assumed  that  the  phase  2  portion  of 
the  ONION-mcthod  algorithm  for  oflhornral  incidence  can 
be  accurately  implemented  in  software  that  treats  the  layers 
of  the  panel  as  fluids.  Once  the  best-fit  properties  have  been 
determined  in  this  manner,  the  phase  3  portion  of  the  algo¬ 
rithm  performs  a  calculation  of  the  reflection  coefficient 
magnitude  as  a  function  of  incidence  angle.  This  calculation 
is  done  by  using  the  best-fit  model  properties  that  have  been 
determined  during  phase  2,  as  well  as  using  assumed  shear 
properties  for  the  layers,  in  a  solid-layer  panel  model.  This 


TABU:  V.  Complex  reflection  coefficients  as  a  function  of  incidence  angle 
for  a  simple  three-layer  sample  Here,  PMM  is  treated  as  a  fluid,  while  steel 
is  treated  as  a  solid.  Frequency  is  20  kHz. 


Incidence 
angle  0 
( degrees ) 

Complex  reflection 
coefficient  (dimensionless) 

Fluid  PMM- 

watcr- 
solid  steel 

0 

0.074,0. 164 

5 

0071.0.182 

in 

0  060.0.2.35 

i? 

0.936.0.320 

20 

0.803.0.424 

25 

0.825.0.544 

30 

0.728.0.668 

35 

0.597.0.785 

40 

0.430.0.883 

45 

0.221.0.946 
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calculation  is  performed  for  an  angular  range  of  0-89  de¬ 
grees  at  increments  of  1°.  (As  previously  mentioned,  the 
0  =  90°  direction  is  omitted  in  order  to  avoid  a  numerical 
singularity  at  tha*  angle  in  the  available  software. )  Comput¬ 
ed  reflection  coefficients  for  incidence  angles  falling  between 
measurement  angles  represent  interpolated  values.  Comput¬ 
ed  reflection  coefficients  for  incidence  angles  greater  than 
the  largest  measured  incidence  angle  represent  extrapolated 
values.  Iu  evaluating  samples  containing  macrovoided  vis¬ 
coelastic  layers,  it  is  only  the  steel  support  plate  that  is  as¬ 
sumed  to  have  a  significant  shear  modulus.  Hence,  negligible 
values  of  the  shear  modulus1"  "  are  used  for  the  macrovoid¬ 
ed  viscoelastic  layers  when  this  final  calculation l:  is  made.  In 
view  of  the  fact  that  this  final  calculation  only  needs  to  be 
performed  once  for  each  measured  frequency,  only  an  insig¬ 
nificant  amount  of  CPU  time  is  required  for  it.  Thus  total 
computation  time  is  not  significantly  increased  by  this  use  of 
the  full  solid-layer  model  in  the  phase  3  calculation. 

If  a  significant  shear-wave  effect,  due  to  the  support 
plate,  arises  at  large  offnormal-incidence  angles,  it  is  as¬ 
sumed  that  this  effect  will  manifest  itself  in  the  final  phase  3 
calculation.  This  idea  can  be  seen  by  examining  Fig.  4(a) 
and  (b).  These  figures  present  graphs  of  the  theoretical  mag¬ 
nitude  of  the  reflection  coefficient  as  a  function  incidence 
angle  for  the  simple  three-layer  panel  discussed  above,  evalu¬ 
ated  here  for  a  20-kHz  interrogating  wave.  Figure  4(a)  pre¬ 
sents  results  based  on  treating  all  three  layers  as  fluids;  Fig¬ 
ure  4(b)  presents  results  based  on  treating  both  the  PMM 
and  steel  layers  as  solids.  The  considerable  influence  of  shear 
waves  in  this  case  can  be  seen  in  the  significant  differences 
between  these  two  graphs  at  the  larger  incidence  angles. 
Note,  however,  that  the  graphs  are  virtually  indistinguish¬ 
able  at  small  incidence  angles.  This  is  an  indication  that  it 
should  be  possible  to  perform  a  least -squares  fit  of  a  fluid- 
layer  model,  as  is  required  in  the  phase  2  portion  of  the  gen¬ 
eralized  ONION  algorithm,  to  data  acquired  from  such  a 
sample,  at  least  for  incidence  angles  that  are  less  than  or 
equal  to  the  incidence  angle  that  corresponds  to  a  rotator 
shaft  angle  of  20°.  After  the  fitting  process  (based  on  the  data 
acquired  at  the  mea  air.  d  angles)  is  complete,  a  curve  closely 
corresponding  to  Fig.  4(b)  (i.e.,  the  solid-layer  curve)  can 
be  generated,  assuming  that  accurate  shear  properties  are 
available  for  each  of  the  layers  of  the  panel.  This  is  what  is 
done  in  the  phase  3  portion  of  the  algorithm. 

III.  EXPERIMENTS 
A.  Simple  three-layer  sample 

In  order  to  provide  a  basis  for  an  initial  investigation  of 
the  effectiveness  of  the  generalized  ONION  method,  experi¬ 
mental  data  were  acquired  from  a  simple  three-layer  sample 
panel.  This  is  the  same  sample  as  was  used  in  the  initial  test¬ 
ing  of  the  normal-incidence  ONION  method.4 '  and  is  the 
same  panel  that  was  considered  theoretically  in  Sec.  II  (C) 
and  (D). 

Data  were  acquired  from  this  sample  at  interrogating- 
wave  frequencies  of  5,  10,  and  20  kHz.  Both  normal  and 
offnormal  incidence  measurements  were  made.  The  rotator- 
shaft  angles  <b  that  were  used1 '  in  these  measurements  are  0°, 
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T«mp:  9.00  D«g  C  Press:  110.00  kPa  Preq:  20.00  kHz 

(•) 


Temp:  9.00  Deg  C  Press:  110.00  kPa  Preq:  20.00  kHz 

(b) 


FIG  4  Theoretical  reflection  coefficient  magnitude  as  a  function  of  inci¬ 
dence  angled  for  the  simple  three*  layer  panel.  Frequency  is  20  kHz.  (a)  All 
panel  layers  treated  as  fluids,  (b)  PMM  and  steel  layers  treated  as  solids. 


5°,  10°,  15°,  and  20°.  A  hydrophone  offset  distance  of  dh  =  5 
cm  from  the  PMM  layer,  and  a  source-to-rotator-shaft  sepa¬ 
ration  distance  off/,  =  200cm,  were  used  in  thisexperiment. 
The  results  of  these  measurements,  and  of  applying  the 


generalized  ONION  method  algorithm  to  the  data,  are  sum¬ 
marized  by  the  graphs  presented  in  Fig.  5.  The  dashed-line 
curves  in  the  graphs  presented  in  these  figures  represent  the 
final  output  of  the  software  based  on  the  phase  3  portion  of 
the  method;  thus  the  dashed-line  curves  were  generated  by 
using  the  best-fit  values  produced  by  the  least-squares  fitting 
process  of  phase  2,  as  well  as  the  known  shear  properties  of 
PMM  and  steel,  in  a  solid-layer  model  calculation.  The  solid 
circles  thai  arc  plotted  on  top  of  these  dashed-line  curves  are 
reference  points  that  are  used  to  indicate  the  incidence  an¬ 
gles  9  that  correspond  to  the  measurem  rnt  angles  d>  at  which 
the  data  were  acquired.  Hence,  portions  of  the  dashed-line 
curves  lying  between  the  solid  circles  represent  interpolated 
reflection  coefficients,  while  portions  of  the  dashed-line 
curves  lying  beyond  the  last  solid  circle  represent  extrapolat¬ 
ed  reflection  coefficients. 

The  effectiveness  of  the  method  in  this  case  can  be  ascer¬ 
tained  by  reference  to  the  solid-line  curves  in  Fig.  5.  These 
curves  were  generated  by  using  the  known  properties  of 
PMM,  water,  and  steel  in  a  theoretical  calculation  based  on 
the  solid-layer  panel  model.  As  can  be  seen  by  comparing  the 
solid-line  and  dashed-line  curves  in  Fig.  5,  the  experimental 
and  theoretical  curves  are  in  good  agreement  at  5  and  10 
kHz,  with  the  exception  of  the  peak  in  the  vicinity  of  0  =  38° 
in  the  5-kHz  curves.  There  is  even  a  reasonable  agreement  in 
the  qualitative  features  of  the  curves  based  on  the  20-kHz 
case,  despite  the  fact  that  this  frequency  might  have  been 
anticipated  to  be  troublesome  for  this  sample.  (Note  again 
Table  IV.)  In  particular,  note  that  the  extrapolation  region 
of  the  dashed-line  curve  presented  in  Fig.  5(c)  (i.e.,  the  re¬ 
gion  beyond  the  last  solid  circle  in  this  figure)  successfully 
predicts  the  “notch"  in  the  vicinity  of  9  =  40°  of  the  20-kHz 
theoretical  curve  (solid  line);  however,  the  reflection  coeffi¬ 
cient  corresponding  to  the  bottom  of  the  notch  is  incorrectly 
predicted  by  approximately  14%.  Quantitative  agreement 
between  the  dashed-line  and  solid-line  curves  of  Fig.  5(c)  is 
poorest  beyond  approximately  9  =  50°,  although  the  quali¬ 
tative  features  are  similar.  Nonetheless,  in  view  of  the  fact 
that  the  data  were  measured  only  up  to  a  rotator-shaft  angle 
of  <f>  =  20°,  the  extrapolation  represents  a  considerable  in¬ 
crease  in  the  amount  of  available  information.  Considering 
the  fact  that  this  panel,  and  particularly  the  20-kHz  frequen¬ 
cy  case,  represents  an  especially  severe  test  of  the  fluid-layer 
panel  model  used  in  the  phase  2  portion  of  the  method,  the 
results  of  this  measurement  can  reasonably  be  described  as 
encouraging. 

B.  Panel  containing  macrovoided  viscoelastic  layers 

Experimental  data  were  next  acquired  from  a  sample 
panel  containing  macrovoided  viscoelastic  layers.  This  is  the 
same  sample  described  in  the  reports  of  the  normal  incider  ;e 
work.4  5 14  A  hydrophone  offset  distance  of  dh  =  5  cm,  and 
source-to-rotator  shaft  separation  distance  of  </,  =  200  cm, 
were  used  in  performing  these  measurements.  The  interro- 
gating-wave  frequencies  were  again  5,  10,  and  20  kHz,  and 
the  hydrostatic  test  pressures  used  were  atmospheric  pres¬ 
sure  and  1380  kPa  (200  psi).  Rotator  shaft  measurement 
angles  <j>  were  0°,  5°,  10°,  15°,  and  20°. 

The  results  of  applying  the  generalized  ONION-meth- 
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MG  5.  Experimental  data  points  (solid  circles)  and  interpolated  and  ex¬ 
trapolated  values  (dashed-line  curves)  compared  to  theory  (solid-line 
curves)  for  the  magnitude  of  the  reflection  coefficient  as  a  function  of  the 
incidence  angle  0  for  the  simple  three-layer  sample,  (a)  5  kHz.  (b)  10kHz. 
and  (0  20  kHz. 


Incidence  Angle  (Degrees) 


od  algorithm  to  these  measurements  are  summarized  in 
Figs.  6  and  7.  Unfortunately,  effective  sound  speeds  and 
losses  of  the  layers  comprised  by  this  sample  are  unknown, 
so  comparison  of  these  results  with  theory  is  not  possible  for 
this  case.  Thus,  in  this  case,  we  examine  instead  the  experi¬ 
mental  and  model  reflected  pulsed  waveforms,  within  the 


appropriate  data  windows,  as  a  function  of  rotator  shaft  an¬ 
gle  6.  Figure  8  presents  these  waveforms  for  the  5-kHz  test 
case  measured  at  the  1380  kPa  (200  psi)  hydrostatic  test 
pressure. ''  As  can  be  seen,  agreement  between  the  model 
and  the  data  is  reasonably  good  at  all  angles.  This  figure  is 
typical  of  the  results  obtained  at  the  other  test  frequencies 
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evaluated  at  this  pressure,  and  is  also  typical  of  the  results 
obtained  for  each  of  the  test  frequencies  evaluated  at  atmo¬ 
spheric  pressure. 

IV.  DISCUSSION 

While  the  experimental  results  presented  for  the  pane! 
containing  simple  homogeneous  layers  are  in  reasonably 


Temp:  9.00  Deg  C  Press:  110.00  kPa  F-eq.  20.00  kHz 
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incidence  Angle  (Degrees) 


FKi  6.  Experimental  data  points  I  solid  circles)  and  interpolated  and  ex¬ 
trapolated  values  (solid-line  curves)  for  the  magnitude  of  the  reflection  co¬ 
efficient  as  a  function  of  incidence  angle  f)  for  the  sample  containing  macro- 
voided  viscoelastic  layers.  Hydrostatic  test  pressure  is  atmospheric,  (a)  5 
kHz.  <b>  10  kHz.  and  <c)  20  kHz 


good  agreement  with  theory,  and  although  the  results  pre¬ 
sented  for  the  panel  containing  macrovoided  viscoelastic 
layers  also  exhibit  reasonable  behavior,  some  further  consid¬ 
erations  are  necessary  to  gain  an  appreciation  of  the  meaning 
of  the  results  obtained  for  these  measurements.  The  ONION 
method  is,  in  essence,  an  extrapolation  procedure;  i.e.,  its 
purpose  is  the  prediction  of  steady-state  results  on  the  basis 
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of  available  transient-wave  data.  In  the  case  of  offnormal 
incidence  panel  measurements,  a  furth-m  extrapolation  to  in¬ 
cidence  angles  that  are  not  experimentally  realizable  (due  to 
the  influence  of  edge  effects)  is  also  done.  In  all  cases,  the 
extrapolations  are  achieved  by  determining  '  litable  model 
parameters  by  performing  a  least-squares  fit  if  a  theoretical 
model  to  experimental  data. 


Tamp:  9.00  Otg  C  Prass:  1380.00  kPo  Freq:  20.00  kHz 

(c) 
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FIG.  7.  Experimental  data  points  (solid  circles)  and  interpolated  and  ex¬ 
trapolated  values  (solid-line  curves)  for  the  magnitude  of  the  reflection  co¬ 
efficient  as  a  function  ->f  incidence  angle  d  for  the  sample  containing  macro- 
voided  viscoelastic  layers.  Hydrostatic  test  pressure  is  1380  kPa  C  X)  psi). 
(a)  5  kHz,  (b)  10  kHz.  and  (c)  20  kHz. 


The  theoretical  model  assumes  effectively  homogeneous 
ers.  Thus  inhomogeneities  within  the  materials  that  form 
the  layers  of  the  sample,  and  discontinuities  represented  by 
the  sample  edges,  are  departures  of  the  experimental  system 
from  that  assumed  by  the  theoretical  model.  In  acquiring 
experimental  data,  s’"ps  are  taken  to  reduce,  as  far  as  possi¬ 
ble,  the  influence  of  experimental  aspects  that  represent  de- 
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DATA  PC'NTS  (EACH  POINT  =  .50  MICROSECONDS) 
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parmres  from  the  model.  One  example  of  this,  previously 
mentioned,  is  the  use  of  narrowing  data  window  widths  at 
increasing  incidence  angles  to  reduce  edge-wave  interfer¬ 
ence 

Each  of  the  macrovoided  viscoelastic  layers  of  test  pan¬ 
els  of  interest  is  usually  not  fabricated  in  the  form  of  a  single, 
continuous  piece  of  material.  Such  layers  are  usually  fabri¬ 
cated  from  a  number  of  subsections  of  material.  (These  sub¬ 
sections  are  called  "tiles.")  These  tiles  typically  have  the 
same  thickness  as  the  layer  to  be  fabricated,  but  have  a 
smaller  cross-sectional  area  than  that  of  the  panel.  (Some¬ 
times  several  tile  thicknesses  are  laminated  together  to  form 
the  desired  layer  thickness.)  Places  where  tiles  are  joined 
together  along  their  edges  in  the  formation  of  a  panel  layer 
are  called  "seams"  Most  panel  designs  use  a  random  tile- 
edge  arrangement  pattern  in  the  fabrication  of  layers.  That 
is,  a  variety  of  tile  cross-sectional  areas  is  used  in  creating  an 
overall  layer  cross  section.  This  is  done  in  order  to  reduce  the 
probability  t hat  the  waves  originating  at  the  tile  seams  will 
constructively  interfere,  (hereby  corrupting  »he  desired  mea¬ 
surement.  However,  neither  the  random  tile  geometry,  nor 
the  use  of  narrowing  data  window  widths,  entirely  elimi¬ 
nates  the  influences  of  unwanted  interfering  waves,  and  the 
presence  of  such  waves  must  be  viewed  as  a  contamination  of 
the  experimental  signal. 

One  consequence  of  the  presence  of  contaminating  edge- 


waves  and  seam  waves  is  that  such  waves  will  result  in  errors 
in  the  determination  of  the  layer  parameters,  i.e.,  sound 
speeds  and  losses.  However,  if  the  mean-squared  error  of  fit 
between  the  model  vvuveioi.n  and  the  experimental  wave¬ 
form  can  be  made  small  (typically,  less  than  10%),  it  is 
assumed  that  such  errors  are  not  significant.  If  this  assump¬ 
tion  is  not  valid,  it  might  very  well  mean  that  any  panel 
measurement  made  on  such  a  sample  is  not  meaningful.  This 
comment  is  based  on  the  fact  that  if,  for  example,  seam  waves 
represent  a  signijicani  component  of  the  observed  wave,  then 
the  experimentally  determined  reflection  coefficient  would 
be  a  function  of  seam  geometry  and  sample  size,  and  this 
would  be  true  regardless  of  the  method  used  to  perform  the 
measurement.  An  extrapolation  of  measurements  obtained 
from  a  panel  having  a  significant  contribution  to  the  ob¬ 
served  siginal  arising  from  a  random  seam  geometry  would 
not  appear  to  be  a  well-posed  problem  and  thus  the  idea  of  a 
panel  measurement  would  not  seem  to  be  meaningful  in  such 
a  ease. 

An  alternate  was  of  v  iewing  the  situation  is  to  imagine 
that  the  observed  wave  represents  the  sum  of:  (i )  a  "seam- 
free  field"  (i.e..  a  field  charaeteriz.able  only  by  the  effective 
medium  properties  of  the  lasers )  and  (ii )  a  "seam  field."  ( In 
the  discussion  that  follows,  we  ignore  the  contribution  from 
panel  edge  waves. )  It  is  expected  that  field  (i )  is  independent 
of  sample  size  and  of  the  details  of  sample  fabrication,  but 
that  field  (ii)  is  dependent  on  each  of  these  factors.  The 
ONION  method  assumes  that  the  influence  of  field  (ii)  is 
negligible,  and  that  field  (i)  can  be  used  to  characterize  the 
panel.  The  extent  to  which  this  assumption  holds  is  tested 
within  the  ONION  method  by  evaluating  the  fit  between  the 
model  and  the  data.  It  is  assumed  that  if  the  "goodness-of- 
lit"  should  ever  happen  to  fall  to  some  unacceptable  level 
(e.g..  worse  than  1?%  )  it  is  due  to  the  influence  of  a  noil- 
negligible  field  of  type  (ii).  In  this  case,  the  extrapolations 
based  on  the  method  would  probably  be  unreliable.  How¬ 
ever,  it  would  seem  that  a  measurement  made  in  such  a  case 
on  an  ostensibly  identical  sample  (i.e.,  one  fabricated  using 
the  same  physical  layer  characteristics  as  the  one  originally 
tested  but  having  a  different  seam  geometry)  would  have  a 
significantly  different  experimental  response  and.  hence,  the 
idea  of  a  panel  measurement  on  such  a  sample  is  probably 
meaningless. 

V.  SUMMARY  AND  CONCLUSIONS 

A  generalization  of  the  ONION-method  algorithm  to 
allow  for  offnormal  incidence  panel  measurements  has  been 
described.  Successful  applications  of  the  technique  to  experi¬ 
mental  measurements  have  also  been  discussed. 

The  generalized  ONION  method  involves  three  phases 
of  calculation.  In  phase  1,  an  ordinary  normal-incidence 
ONION  calculation  determines  initial  model  parameters;  in 
phase  2,  a  simultaneous  least-squares  fit  of  an  ofTnormal  inci¬ 
dence  fluid-layer  panel  model  to  experimental  data  acquired 
over  a  range  of  incidence  angles  is  performed;  in  phase  3,  the 
best-fit  properties  deduced  by  phase  2,  and  a  priori  shear 
properties  for  the  panel  layers,  are  used  in  a  solid-layer  panel 
model  calculation  to  obtain  interpolated  and  extrapolated 
reflection  coefficients  as  a  function  of  incidence  angle. 
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The  experimental  measurements  that  were  described 
here  involved  a  simple  three-layer  sample  as  well  as  a  sample 
containing  macrovoided  viscoelastic  layers.  The  simple 
three-layer  sample  represents  a  severe  test  of  the  fluid-layer 
panel  model  used  in  phase  2,  due  to  the  significant  shear 
modulus  and  thickness  of  the  TMivi  iayer.  Agreement 
between  measurements  obtained  from  this  sample  and  theo¬ 
retical  calculations  was  found  to  be  reasonably  good.  The 
results  obtained  from  the  sample  containing  macrovoided 
viscoelastic  layers  also  exhibit  reasonable  behavior. 

The  generalized  ONION  method  assumes  that  the  lay¬ 
ers  of  the  sample  panel  can  be  characterized  using  effective 
medium  properties.  If  actual  sample  inhomogeneities  such 
as  seams  and  edges  produce  substantial  contributions  to  the 
observed  field,  the  extrapolations  will  probably  be  unrelia¬ 
ble;  however,  in  this  case,  a  panel  measurement  is  probably 
not  meaningful. 

The  present  paper  represents  an  initial  investigation  into 
the  generalization  of  the  ONION  method  to  offnormal  inci¬ 
dence.  Future  work  will  seek  to  establish  well-defined  re¬ 
gions  of  applicability  of  the  generalized  method,  and  will 
further  evaluate  the  possibility  of  using  an  effective  medium 
theory  to  more  accurately  account  for  the  influence  of  shear 
waves. 
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Am.  63.  1694-1698  (1978). 

"Instead  of  simply  assuming  a  negligible  shear  modulus,  it  may  be  possible 
to  use  an  effective  shear  modulus  for  the  macrovoided  v  iscoelastic  layers 
based  on  one  of  the  many  available  effective  medium  theories.  For  an  ex¬ 
ample  of  one  of  these,  see  the  summary  provided  in  G.  C.  Gaunaurd  and 
W  Wcriman.  "Comparison  of  effective  medium  theories  for  inhomogen¬ 
eous  continua."  J  Acoust  Soc.  Am  85.  541-554  (  1989).  However,  in 
v  lew  of  the  fact  that  effective  medium  theories  are  somewhat  controversial 
|  see  the  article,  which  accompanies  the  one  mentioned  above,  by  L.  W 
Anson  and  R.  C.  Olivers.  "Ultrasonic  propagation  in  suspensions.  A 
comparison  of  a  multiple  scattering  and  an  effective  medium  approach." 
J.  Acoust.  Soc.  Am.  85.  535-540  ( 1989) ).  and  ill  view  the  fact  that  vali¬ 
dating  data  are  scarce.  the  possible  use  of  an  effective  medium  theory  is  left 
as  an  open  issue  here. 

'  Since  the  PMM  layer  of  the  simple  three-layer  sample  has  significant 
shear  properties,  the  actual  value  of  the  shear  modulus  of  PMM  must  be 
used  in  implementing  phase  3  for  measurements  made  on  this  sample. 
However,  this  presents  no  difficulty  since  the  shear  modulus  of  PMM  is 
known. 

'  As  discussed  in  Sec.  II.  these  rotator  shaft  angles  i  must  be  corrected 
using  Eq.  ( 1 )  to  obtain  the  incidence  angles  0 
"This  panel  is  composed  of  three  macrovoided  viscoelastic  layers  and  a 
fourth  layer  that  isa().95-cm  ( xin.  I  steel  backing  plate.  The  first  layer  has 
5-cm  thickness,  the  second  layer  4-em  thickness,  and  the  ihird  layer  3-ctn 
thickness.  The  acoustical  properties  of  the  lay  ers  are  unknown. 

"The  endpoints  of  the  data  windows  depicted  in  Fig.  8  were  chosen  to  be 
consistent  with  Eq  (4).  However,  the  start  points  of  these  data  windows 
avoid  the  first  300  data  points  permitted  by  Eq.  (4).  This  is  done  to  ac¬ 
count  for  the  noiicausal  nature  of  the  ONION  model.  (See  Ref.  4. )  Note 
that  the  number  of  data  points  given  for  each  of  the  windows  in  Fig.  8 
corresponds  to  a  data  measured  rale  of  2  MHz.  which  differs  from  the  4 
MHz  data  measurement  rate  used  in  the  example  given  in  connection  with 
Eq,  (4). 


1427 


J.  Acoust.  Soc.  Am.,  Vol.  87,  No.  4,  April  1990 


Jean  C.  Piquette:  Offnormal  reflection  coefficient 


1427 


